Recent Theoretical Results for Nonequilibrium 
Deposition of Submicron Particles 



Vladimir Privman 



Department of Physics, Clarkson University, Potsdam, New York 13699-5820, USA 



ABSTRACT 

Selected theoretical developments in modeling of deposition of sub- 
micrometer size (submicron) particles on solid surfaces, with and without 
surface diffusion, of interest in colloid, polymer, and certain biological sys- 
tems, are surveyed. We review deposition processes involving extended 
objects, with jamming and its interplay with in-surface diffusion yield- 
ing interesting dynamics of approach to the large-time state. Mean-field 
and low-density approximation schemes can be used in many instances 
for short and intermediate times, in large enough dimensions, and for 
particle sizes larger than few lattice units. Random sequential adsorp- 
tion models are appropriate for higher particle densities (larger times). 
Added diffusion allows formation of denser deposits and leads to power- 
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law large-time behavior which, in one dimension (linear substrate, such as 
DNA), was related to diffusion- limited reactions, while in two dimensions 
(planar substrate) , was associated with evolution of the domain- wall and 
defect network, reminiscent of equilibrium ordering processes. 



This is a review article, to appear in The Journal of Adhesion (2000). 
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1. Introduction 

1.1. Surface Deposition of Submicron Particles 

Surface deposition of submicron particles is of immense practical im- 
portance [1-4]. Typically, particles of this size, colloid, protein or other 
biological objects, are suspended in solution, without sedimentation due 
to gravity. In order to maintain the suspension stable, one has to prevent 
aggregation (coagulation) that results in larger flocks for which gravity 
pull is more profound. Stabilization by particle-particle electrostatic re- 
pulsion or by steric effects, etc., is usually effective for a sufficiently dilute 
suspension. But this means that even if a well-defined suspension of well- 
characterized particles is available, it cannot be always easily observed 
experimentally in the bulk for a wide range of particle interactions. For 
those interaction parameters for which the system is unstable with respect 
to coagulation, the time of observation will be limited by the coagulation 
process which can be quite fast. 

One can form a dense deposit slowly, if desired, on a surface. Indeed, 
particles can be deposited by diffusion, or more realistically by convective 
diffusion [5] from a flowing suspension, on collector surfaces. The suspen- 
sion itself need not be dense even though the on-surface deposit might be 
quite dense, depending of the particle-particle and particle-surface inter- 
actions. Dilution of suspension generally prolongs an experiment aimed 
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at reaching a certain surface coverage. Thus, surface deposition has been 
well established as an important tool to probe interactions of matter ob- 
jects on the submicron scale [1-4]. 

1.2. Particle Jamming and Screening at Surfaces 

Figure 1 illustrates possible configurations of particles at a surface. 
From left to right, we show particles deposited on the surface of a collec- 
tor, then particles deposited on top of other particles. The latter is possi- 
ble only in the absence of significant particle-particle repulsion. The two 
situations are termed monolayer and multilayer deposition even though 
the notion of a layer beyond the one exactly at the surface is only ap- 
proximate. We next show two effects that play important role in surface 
growth. The first is jamming: a particle marked by an open circle cannot 
fit in the lowest layer at the surface. A more realistic two-dimensional 
(2D) configuration is shown in the inset. 

The second effect is screening: surface position marked by the open 
circle is not reachable. Typically, in colloid deposition monolayer or few- 
layer deposits are formed and the dominant effect is jamming, as will 
be discussed later. Screening plays dominant role in deposition of mul- 
tiple layers and, together with the transport mechanism, determines the 
morphology of the growing surface. In addition, the configuration on the 
surface depends on the transport mechanism of the particles to it and 
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on the particle motion on the surface, as well as possible detachment. 
Particle motion is typically negligible for colloidal particles but may be 
significant for proteins. 

1.3. Role of Dimensionality and Relation to Other Systems 

An important feature of surface deposition is that for all practical 
purposes it is essentially a 2D problem. As a result, any mean-field, rate- 
equation, effective-field, etc., approaches which are usually all related 
in that they ignore long-range correlations and fluctuation effects, may 
not be applicable. Indeed, it is known that as the dimensionality of a 
many-body interacting system decreases, fluctuations play a larger role. 
Dynamics of important physical, chemical, and biological processes [6-7] 
provides examples of strongly fluctuating systems in low dimensions, D = 
1 or 2. These processes include surface adsorption on planar substrates 
or on large collectors. The surface of the latter is semi- two-dimensional 
owing to their large size as compared to the size of the deposited particles. 

The classical chemical reaction-diffusion kinetics corresponds to D = 
3. However, heterogeneous catalysis generated interest in D = 2. For 
both deposition and reactions, some experimental results exist even in 
D = 1 (see later). Finally, kinetics of ordering and phase separation, 
largely amenable to experimental probe in D = 3 and 2, attracted much 
recent theoretical effort in D = 1, 2. 
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Models m D = 1, and sometimes in D = 2, allow derivation of an- 
alytical results. Furthermore, it turns out that all three types of model: 
deposition-relaxation, reaction-diffusion, phase separation, are interre- 
lated in many, but not all, of their properties. This observation is by no 
means obvious. It is model-dependent and can be firmly established [6-7] 
only in low dimensions, mostly in D = 1. 

Such low-dimensional nonequilibrium models pose several interesting 
challenges theoretically and numerically. While many exact, asymptotic, 
and numerical results are already available in the literature [6-7], this 
field presently provides examples of properties which lack theoretical ex- 
planation even in ID. Numerical simulations are challenging and require 
large scale computational effort already for ID models. For more experi- 
mentally relevant 2D cases, where analytical results are scarce, difficulty 
in numerical simulations has been the limiting factor in understanding of 
many open problems. 

1.4- Outline of This Review 

The purpose of this article is to provide an introduction to the field 
of nonequilibrium surface deposition models of extended particles. By 
"extended" we mean that the main particle-particle interaction effect 
will be jamming, i.e., mutual exclusion. No comprehensive survey of the 
literature is attempted. Relation of deposition to other low-dimensional 



-6- 



models mentioned earlier will be only referred to in detail in few cases. 
The specific models and examples selected for a more detailed exposi- 
tion, i.e., models of deposition with diffusional relaxation, were biased by 
author's own work. 

The outline of the review is as follows. The rest of this introductory 
section is devoted to defining the specific topics of surface deposition 
to be surveyed. Section 2 describes the simplest models of random se- 
quential adsorption. Section 3 is devoted to deposition with relaxation, 
with general remarks followed by definition of the simplest, ID models 
of diffusional relaxation for which we present a more detailed description 
of various theoretical results. Multilayer deposition is also commented 
on in Section 3. More numerically-based 2D results for deposition with 
diffusional relaxation are surveyed in Section 4. Section 5 presents brief 
concluding remarks. 

Surface deposition is a vast field of study. Our emphasis here will 
be on those deposition processes where the particles are "large" as com- 
pared to the underlying atomic and morphological structure of the sub- 
strate and as compared to the range of the particle-particle and particle- 
substrate interactions. Thus, colloids, for instance, involve particles of 
submicron to several micron size. We note that l/j,m= 10000 A, whereas 
atomic dimensions are of order lA, while the range over which particle- 
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surface and particle-particle interactions are significant as compared to 
fcT, is typically of order lOOA or less. Extensive theoretical study of such 
systems is relatively recent and it has been motivated by experiments 
where submicron-size colloid, polymer, and protein "particles" were the 
deposited objects [1-4,8-18]. 

Perhaps the simplest and the most studied model with particle ex- 
clusion is Random Sequential Adsorption (RSA). The RSA model, to be 
described in detail in Section 2, assumes that particle transport (incom- 
ing flux) onto the surface results in a uniform deposition attempt rate R 
per unit time and area. In the simplest formulation, one assumes that 
only monolayer deposition is allowed. Within this monolayer deposit, 
each new arriving particle must either fit in an empty area allowed by 
the hard-core exclusion interaction with the particles deposited earlier, 
or the deposition attempt is rejected. 

The basic RSA model will be described shortly, in Section 2. Recent 
work has been focused on its extensions to allow for particle relaxation 
by diffusion, see Sections 3 and 4, to include detachment processes, and 
to allow multilayer formation. The latter two extensions will be briefly 
surveyed in Section 3. Several other extensions will not be discussed [1-4]. 
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2. Random Sequential Adsorption 

2.1. The RSA Model 

The irreversible Random Sequential Adsorption (RSA) process [19- 
20] models experiments of submicron particle deposition by assuming a 
planar 2D substrate and, in the simplest case, continuum (off-lattice) 
deposition of spherical particles. However, other RSA models have also 
received attention. In 2Z), noncircular cross-section shapes as well as 
various lattice-deposition models were considered [19-20]. Several experi- 
ments on polymers and attachment of fluorescent units on DNA molecules 
[18] (the latter is usually accompanied by motion of these units on the 
DNA and detachment) suggest consideration of the lattice-substrate RSA 
processes in ID. RSA processes have also found applications in traffic 
problems and certain other fields. Our presentation in this section aims 
at defining some RSA models and outlining characteristic features of their 
dynamics. 

Figure 2 illustrates the simplest possible monolayer lattice RSA model: 
irreversible deposition of dimers on the linear lattice. An arriving dimer 
will be deposited if the underlying pair of lattice sites are both empty. 
Otherwise, it is discarded, which is shown schematically by the two dimers 
above the surface layer. Their deposition on the surface is not possible 
unless detachment and/or motion of monomers or whole dimers clear the 
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appropriate landing sites. 

Let us consider the irreversible RSA without detachment or diffusion. 
The substrate is usually assumed to be empty initially, at t = 0. In the 
course of time t, the coverage, p{t), increases and builds up to order 1 
on the time scales of order {RV)~^, where R was defined earlier as the 
deposition attempt rate per unit time and area of the surface, while V 
is the particle D-dimensional "volume." For deposition of spheres on a 
planar surface, V is actually the cross-sectional area. 

At large times the coverage approaches the jammed-state value where 
only gaps smaller than the particle size were left in the monolayer. The 
resulting state is less dense than the fully ordered close-packed coverage. 
For the D = 1 deposition shown in Figure 2 the fully ordered state would 
have p = 1 . The variation of the RSA coverage is illustrated by the lower 
curve in Figure 3. 

At early times the monolayer deposit is not dense and the deposition 
events are largely uncorrelated. In this regime, mean-field like low-density 
approximation schemes are useful [21-23]. Deposition of k-mei particles 
on the linear lattice in ID was in fact solved exactly for all times [24]. 
In D = 2, extensive numerical studies were reported [23,25-36] of the 
variation of coverage with time and large-time asymptotic behavior which 
will be discussed shortly. Some exact results [24] for correlation properties 
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are available in ID. Numerical results [27] for correlation properties have 
been obtained in 2D. 

2.2. The Large- Time Behavior in RSA 

The large-time deposit has several characteristic properties. For lat- 
tice models, the approach to the jammed-state coverage is exponential 
[36-38] . This was shown to follow from the property that the final stages 
of deposition are in few sparse, well separated surviving landing sites. 
Estimates of decrease in their density at late stages suggest that 



p(oo) - p{t) ~ exp {-m^t) , (1) 

where I is the lattice spacing and D is the dimensionality of the substrate. 
The coefficient in Eq. (1) is of order /V ii the coverage is defined as the 
fraction of lattice units covered, i.e., the dimensionless fraction of area 
covered, also termed the coverage fraction, so that coverage as density 
of particles per unit volume would be V~^p. The detailed behavior de- 
pends of the size and shape of the depositing particles as compared to 
the underlying lattice unit cells. 

However, for continuum off-lattice deposition, formally obtained as 
the limit ^ — ^ 0, the approach to the jamming coverage is power-law. 
This interesting behavior [37-38] is due to the fact that for large times the 



remaining voids accessible to particle deposition can be of sizes arbitrarily 
close to those of the depositing particles. Such voids are thus reached 
with very low probability by the depositing particles, the flux of which 
is uniformly distributed. The resulting power-law behavior depends on 
the dimensionality and particle shape. For instance, for Z)-dimensional 
cubes of volume V, 

while for spherical particles. 
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p(oo) - Pit) ~ (RVt)-^/" . (3) 

For D > 1, the expressions Eqs. (2-3), and similar relations for other 
particle shapes, are actually empirical asymptotic laws which have been 
verified, mostly for D = 2, by extensive numerical simulations [4,25- 
36]. The most studied 2D geometries are circles (corresponding to the 
deposition of spheres on a planar substrate) and squares. The jamming 
coverages are 



Psquares(oo) ~ 0.5620 and Pcircles(oo) C:^ 0.544 to 0.550 , (4) 
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much lower than the close-packing values, 1 and ~ 0.907, respec- 
tively. For square particles, the crossover to continuum in the limit 
k ^ oo and i ^ 0, with fixed V^^^ = ki in deposition oikxkx. . .xk lat- 
tice squares, has been investigated in some detail [36], both analytically 
(in any D) and numerically (in 2D). 

The correlations in the large-time jammed state are different from 
those of the equilibrium random gas of particles with density near p{oo). 
In fact, the two-particle correlations in continuum deposition develop a 
weak singularity at contact, and correlations generally reflect the full 
irreversibility of the RSA process [24,27,38]. 
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3. Deposition with Relaxation 

3.1. Detachment and Diffusional Relaxation 

Monolayer deposits may relax, i.e., explore more configurations, by 
particle motion on the surface, by their detachment, as well as by mo- 
tion and detachment of the constituent monomers or recombined units. 
In fact, detachment has been experimentally observed in deposition of 
colloid particles which were otherwise quite immobile on the surface [39] . 
Theoretical interpretation of colloid particle detachment data has proved 
difficult, however, because binding to the substrate once deposited, can 
be different for different particles, whereas the transport to the substrate, 
i.e., the flux of the arriving particles in the deposition part of the process, 
typically by convective diffusion, is more uniform. Detachment also plays 
role in deposition on DNA molecules [18]. 

Recently, more theoretically motivated studies of the detachment re- 
laxation processes, in some instances with surface diffusion allowed as 
well, have lead to interesting model studies [40-46]. These investigations 
did not always assume detachment of the original units. Models involving 
monomer recombination prior to detachment, of /c-mers in D = 1, have 
been mapped onto certain spin models and symmetry relations were iden- 
tified which allowed derivation of several exact and asymptotic results on 
the correlations and other properties [40-46]. We note that deposition 
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and detachment combine to drive the dynamics into a steady state, rather 
than jammed state as in ordinary RSA. These studies have been largely 
limited thus far to ID models. 

We now turn to particle motion on the surface, in a monolayer de- 
posit, which was experimentally observed in deposition of proteins [17] 
and also in deposition on DNA molecules [18,47]. From now on, we con- 
sider diffusional relaxation, i.e., random hopping on the surface in the 
lattice case. The dimer deposition in ID, for instance, is shown in Fig- 
ure 2. Hopping of dimer particles one site to the left or to the right is 
allowed only if the target site is not occupied. Such hopping can open a 
two-site gap to allow additional deposition. Thus, diffusional relaxation 
lets the deposition process to reach denser, in fact, close-packed config- 
urations. Initially, for short times, when the empty area is plentiful, the 
effect of the in-surface particle motion will be small. However, for large 
times, the density will exceed that of the RSA process, as illustrated by 
the upper curve in Figure 3. 

It is important to emphasize that deposition and diffusion are two 
independent processes going on at the same time. External particles 
arrive at the surface with a fixed rate per unit area. Those finding open 
landing sites are deposited; others are discarded. At the same time, 
internal particles, those already on the surface, attempt, with some rate. 
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to hop to a nearby site. They actually move only if the target site is 
available. 

3.2. One- Dimensional Models 

Further investigation of this effect is much simpler in ID than in 2D. 
Let us therefore consider the ID case first, postponing the discussion of 
2D models to the next section. Specifically, consider deposition of k-meis 
of fixed length V. By keeping the length fixed, we can also naturally 
consider the continuum limit of no lattice by having the lattice spacing 
vanish as — ^ oo. This limit corresponds to continuum deposition if we 
take the underlying lattice spacing £ = V/k. Since the deposition attempt 
rate R was defined per unit area (unit length here), it has no significant 
/c-dependence. However, the added diffusional hopping of /c-mers on the 
ID lattice, with the attempt rate to be denoted by H, and hard-core 
or similar particle interaction, must be fc-dependent. Indeed, we consider 
each deposited k-mer particle as randomly and independently attempting 
to move one lattice spacing to the left or to the right with the rate H/ 2 per 
unit time. Particles cannot run over each other so some sort of hard-core 
interaction must be assumed, i.e., in a dense state most hopping attempts 
will fail. However, if left alone, each particle would move diffusively for 
large time scales. In order to have the resulting diffusion constant V 
finite in the continuum limit — ^ oo, we must assume that 
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(5) 



which is only vahd in ID. 

Each successful hopping of a particle results in motion of one empty 
lattice site. It is useful to reconsider the dynamics of particle hopping 
in terms of the dynamics of this rearrangement of empty area fragments 
[48-50] . Indeed, if several of these empty sites are combined to form large 
enough voids, deposition attempts can succeed in regions of particle den- 
sity which would be jammed in the ordinary RSA. In terms of these new 
"diffuser particles" which are the empty lattice sites of the deposition 
problem, the process is in fact that of reaction-diffusion. Indeed, k reac- 
tants (empty sites) must be brought together by diffusional hopping in 
order to have finite probability of their annihilation, i.e., disappearance 
of a group of consecutive nearest-neighbor empty sites due to successful 
deposition. Of course, the fc-group can also be broken apart due to dif- 
fusion. Therefore, the fc-reactant annihilation is not instantaneous in the 
reaction nomenclature. Such A;-particle reactions are of interest on their 
own [51-57]. 

3.3. Beyond the Mean-Field Approximation 

The simplest mean-field rate equation for annihilation of k reactants 
describes the time dependence of the coverage, p(t), in terms of the reac- 



tant density 1 — p, i.e., the density of the empty spaces, 



dp 

lit 



r(l-p) 



k 



(6) 



where V is the effective rate constant. Note that we assume that the 
close-packing dimensional coverage is 1 in ID. There are two problems 
with this approximation. Firstly, it turns out that for = 2 the mean- 
field approach breaks down. Diffusive-fluctuation arguments for non- 
mean-field behavior have been advanced for several chemical reactions 
[51,53,59-59]. In ID, several exact calculations support this conclusion 
[60-66]. The asymptotic large-time behavior turns out to be 



rather than the mean-field prediction ^ 1/t. The coefficient in Eq. (7) is 
expected to be universal, when expressed in an appropriate dimensionless 
form by introducing single-reactant diffusion constant. 

The power law Eq. (7) was confirmed by extensive numerical sim- 
ulations of dimer deposition [67] and by exact solution [68] for one par- 
ticular value of H for a model with dimer dissociation. The latter work 
also yielded some exact results for correlations. Specifically, while the 
connected particle-particle correlations spread diffusively in space, their 
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p^l/Vi {k = 2,D = l) 



(7) 



decay it time is nondiffusive [68]. Series expansion studies of models of 
dimer deposition with diffusional hopping of the whole dimers or their 
dissociation into hopping monomers, has confirmed the expected asymp- 
totic behavior and also provided estimates of the coverage as a function 
of time [69]. 

The case A; = 3 is marginal with the mean-field power law modified 
by logarithmic terms. The latter were not observed in Monte Carlo stud- 
ies of deposition [49]. However, extensive results are available directly 
for three-body reactions [53-56], including verification of the logarithmic 
corrections to the mean-field behavior [54-56]. 

3.4- Continuum Limit of Off-Lattice Deposition 

The second problem with the mean-field rate equation is identified 
when one attempts to use it in the continuum limit corresponding to 
off-lattice deposition, i.e., for k ^ oo. Note that Eq. (6) has no regular 
limit as — ^ oo. The mean-field approach is essentially the fast diffusion 
approximation assuming that diffusional relaxation is efficient enough to 
equilibrate nonuniform density fluctuations on the time scales fast as com- 
pared to the time scales of the deposition events. Thus, the mean-field 
results are formulated in terms of the uniform properties, such as the den- 
sity. It turns out, however, that the simplest, k^^-powei of the reactant 
density form Eq. (6) is only appropriate for times t >> e'^'^ / (RV). 
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This conclusion was reached [48] by assuming the fast-diffusion, ran- 
domized hard-core reactant system form of the inter-reactant distribution 
function in ID. This approach, not detailed here, allows estimation of 
the limits of validity of the mean-field results and it correctly suggests 
mean-field validity for /c = 4, 5, . . ., with logarithmic corrections for k = 3 
and complete breakdown of the mean-field assumptions for = 2. This 
detailed analysis yields the modified mean-field relation 

dp ^RV(l-pf 

Tt - (1-p + fc-V) (D-l) , (8) 

where 7 is some effective dimensionless rate constant. This new expres- 
sion applies uniformly as /c ^ 00. Thus, the continuum deposition is also 
asymptotically mean-field, with the essentially-singular rate equation 



^ = 7(1 - P) exp[-p/(l - p)] {k = oc,D = l) . (9) 

The approach to the full, saturation coverage for large times is extremely 
slow. 



Similar predictions were also derived for /c-particle chemical reactions 
[53]. 
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3. 5. Comments on Multilayer Deposition 

When particles are allowed to attach also on top of each other, with 
possibly some rearrangement processes allowed as well, multilayer de- 
posits will be formed. It is important to note that the large-layer struc- 
ture of the deposit and fluctuation properties of the growing surface will 
be determined by the transport mechanism of particles to the surface and 
by the allowed relaxations (rearrangements). Indeed, these two charac- 
teristics determine the screening properties of the multilayer formation 
process which in turn shape the deposit morphology, which can range 
from fractal to dense, and the roughening of the growing deposit sur- 
face. There is a large body of research studying such growth, with recent 
emphasis on the growing surface fluctuation properties. 

However, the feature characteristic of the RSA process, i.e., the ex- 
clusion due to particle size, plays no role in determining the universal, 
large-scale properties of thick deposits and their surfaces. Indeed, the 
RSA- like jamming will be only important for detailed morphology of the 
first few layers in a multilayer deposit. However, it turns out that RSA- 
like approaches (with relaxation) can be useful in modeling granular com- 
paction [70]. 

In view of the above remarks, multilayer deposition models involving 
jamming effects were relatively less studied. They can be divided into two 
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groups. Firstly, structure of the deposit in the first few layers is of interest 
[71-73] because they retain memory of the surface. Variation of density 
and other correlation properties away from the wall has structure on the 
length scale of particle size. These typically oscillatory features decay 
away with the distance from the wall. Numerical Monte Carlo simulation 
aspects of continuum multilayer deposition (ballistic deposition of SD 
balls) were reviewed in [73]. Secondly, few-layer deposition processes 
have been of interest in some experimental systems. Mean-field theories 
of multilayer deposition with particle size and interactions accounted for 
were formulated [74] and used to fit such data [15-16,75-76]. 
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4. Two-Dimensional Deposition with Diffusional Relaxation 

4.1. Combined Effects of Jamming and Diffusion 

We now turn to the 2D case of deposition of extended objects on 
planar surfaces, accompanied by diffusional relaxation, assuming mono- 
layer deposits. We note that the available theoretical results are limited 
to few studies [34,77-79]. They indicate a rich pattern of new effects as 
compared to ID. In fact, there exists extensive literature [81] on de- 
position with diffusional relaxation in other models, in particular those 
where the jamming effect is not present or plays no significant role. These 
include deposition of monomer particles, usually of atomic dimensions, 
which align with the underlying lattice without jamming, as well as mod- 
els where many layers are formed (mentioned in the preceding section). 

The 2D deposition with relaxation of extended objects is of interest 
in certain experimental systems where the depositing objects are proteins 
[17]. Here we focus on the combined effect of jamming and diffusion, and 
emphasize dynamics at large times. For early stages of the deposition 
process, low-density approximation schemes can be used. One such ap- 
plication was reported [34] for continuum deposition of circles on a plane. 

In order to identify new features characteristic of 2D, let us consider 
deposition of 2 x 2 squares on the square lattice. The particles are exactly 
aligned with the 2x2 lattice sites as shown in Figure 4. Furthermore, 
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we assume that the diffusional hopping is along the lattice directions ±x 
and ±y, one lattice spacing at a time. In this model dense configurations 
involve domains of four phases as shown in Figure 4. As a result, im- 
mobile fragments of empty area can exist. Each such single-site vacancy 
(Figure 4) serves as a meeting point of four domain walls. 

Here by "immobile" we mean that the vacancy cannot move due to 
local motion of the surrounding particles. For it to move, a larger empty- 
area fragment must first arrive, along one of the domain walls. One such 
larger empty void is shown in Figure 4. Note that it serves as a kink 
in the domain wall. Existence of locally immobile ("frozen") vacancies 
suggests possible frozen glassy behavior with extremely slow relaxation, 
at least locally. The full characterization of the dynamics of this model 
requires further study. The first numerical results [77] do provide some 
answers which will be reviewed shortly. 

4-2. Ordering by Shortening of Domain Walls 

We first consider a simpler model depicted in Figure 5. In this model 
[78-79] the extended particles are squares of size \/2 x y/2. They are 
rotated 45° with respect to the underlying square lattice. Their diffusion, 
however, is along the vertical and horizontal lattice axes, by hopping 
one lattice spacing at a time. The equilibrium variant of this model 
(without deposition, with fixed particle density) is the well-studied hard- 
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square model [82] which, at large densities, phase separates into two 
distinct phases. These two phases also play role in the late stages of 
RSA with diffusion. Indeed, at large densities the empty area is stored in 
domain walls separating ordered regions. One such domain wall is shown 
in Figure 5. Snapshots of actual Monte Carlo simulation results can be 
found in Refs. 78-79. 

Figure 5 illustrates the process of ordering which essentially amounts 
to shortening of domain walls. In Figure 5, the domain wall gets shorter 
after the shaded particles diffusively rearrange to open up a deposition 
slot which can be covered by an arriving particle. Numerical simulations 
[78-79] find behavior reminiscent of the low-temperature equilibrium or- 
dering processes [83-85] driven by diffusive evolution of the domain-wall 
structure. For instance, the remaining uncovered area vanishes according 
to 



1 - p{t) ~ . (11) 

This quantity, however, also measures the length of domain walls in the 
system (at large times). Thus, disregarding finite-size effects and assum- 
ing that the domain walls are not too convoluted (as confirmed by numer- 
ical simulations), we conclude that the power law Eq. (11) corresponds to 
typical domain sizes growing as ~ reminiscent of the equilibrium or- 



dering processes of systems with nonconserved order parameter dynamics 
[83-85]. 

4.3. Numerical Results for Models with Frozen Vacancies 

We now turn again to the 2x2 model of Figure 4. The equihbrium 
variant of this model corresponds to hard-squares with both nearest and 
next-nearest neighbor exclusion [82,86-87]. It has been studied in lesser 
detail than the two-phase hard-square model described in the preceding 
paragraphs. In fact, the equilibrium phase transition has not been fully 
classified (while it was Ising for the simpler model). The ordering at low 
temperatures and high densities was studied [86]. However, many features 
noted, for instance large entropy of the ordered arrangements, require 
further investigation. The dynamical variant (RSA with diffusion) of 
this model was studied numerically [77] . The configuration of the single- 
site frozen (locally immobile) vacancies and the associated network of 
domain walls turn out to be boundary-condition sensitive. For periodic 
boundary conditions the density freezes at values 1 — p ~ L~^, where L 
is the linear system size. 

Preliminary indications were found [77] that the domain size and 
shape distributions in such a frozen state are nontrivial. Extrapolation 
L — ^ 00 indicates that the power law behavior similar to Eq. (11) is 
nondiffusive: the exponent 1/2 is replaced by ~ 0.57. However, the 
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density of the smallest mobile vacancies, i.e., dimer kinks in domain walls, 
one of which is illustrated in Figure 4, does decrease diffusively. Further 
studies are needed to fully clarify the ordering process associated with 
the approach to the full coverage as t — ^ oo and L — ^ oo in this model. 

Even more complicated behaviors are possible when the deposit- 
ing objects are not symmetric and can have several orientations as they 
reach the substrate. In addition to translational diffusion (hopping), one 
has to consider possible rotational motion. The square-lattice deposition 
of dimers, with hopping processes including one-lattice-spacing motion 
along the dimer axis and 90° rotations about a constituent monomer, 
was studied [80] . The dimers were allowed to deposit vertically and hor- 
izontally. In this case, the full close-packed coverage is not achieved at 
all because the frozen vacancy sites can be embedded in, and move by 
diffusion in, extended structures of different topologies. These structures 
are probably less efficiently demolished by the motion of mobile vacan- 
cies than the elimination of localized frozen vacancies in the model of 
Figure 4. 
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5. Conclusion 

In summary, we reviewed theoretical developments in the description 
of deposition processes of extended objects, with jamming and diffusional 
relaxation. While significant progress has been achieved in ID, the 2D 
systems require further study. Most of these investigations will involve 
large-scale numerical simulations. 

Other research directions that require further work include multi- 
layer deposition and particle detachment, especially the theoretical de- 
scription of the latter, including the description of the distribution of 
values/shapes of the primary minimum in the particle-surface interaction 
potential. This would allow to advance beyond the present theoretical 
trend of studying deposition as mainly the process of particle transport 
to the surface, with little or no role played by the details of the actual 
particle-surface and particle-particle double-layer and other interactions. 
Ultimately, we would like to interrelate the present deposition studies and 
approaches in study of adhesion [4], of typically larger particles of sizes 
up to several microns, at surfaces. 
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Figure Captions 

Figure 1: Possible configurations of particles at surfaces. From left to 
right, A — particles deposited directly on the collector; B — parti- 
cles deposited on top of other particles. We next show an example 
of jamming, C — a particle marked by an open circle cannot fit in 
the lowest layer at the surface. A top view of a more realistic two- 
dimensional (2D) surface configuration is shown in the inset. The 
rightmost example, E, illustrates screening: surface position marked 
by the open circle is not reachable. 

Figure 2: Deposition of dimers on the ID lattice. Only one of the three 
hatched dimers can deposit on the surface, which then becomes fully 
jammed in the interval shown. 

Figure 3: Schematic variation of the coverage p{t) with time for deposition 
without (lower curve) and with (upper curve) diffusional or other 
relaxation. The "ordered" density corresponds to close packing. 

Figure 4: Fragment of a deposit configuration in the deposition of 2 x 2 
squares. Illustrated are one single-site frozen vacancy at which four 
domain walls meet (indicated by arrows), and one dimer vacancy 
which causes a kink in one of the domain walls. 

Figure 5: Illustration of deposition of x a/2 particles on the square 
lattice. Diffusional motion during time interval from ti to t2 can 



rearrange the empty area "stored" in the domain wall to open up 
a new landing site for deposition. This is illustrated by the shaded 
particles. 
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